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Abstract. An adiabatic, non-radial pulsation study of a 0.563 M©DA white dwarf model is presented on the basis 
of new evolutionary calculations performed in a self-consistent way with the predictions of time dependent element 
diffusion, nuclear burning and the history of the white dwarf progenitor. Emphasis is placed on the role played by 
the internal chemical stratification of these new models in the behaviour of the eigenmodes, and the expectations 
for the full g-spectrum of periods. The implications for the mode trapping properties are discussed at length. 
In this regard, we find that, for high periods, the viability of mode trapping as a mode selection mechanism is 
markedly weaker for our models, as compared with the situation in which the hydrogen-helium transition region 
is treated assuming equilibrium diffusion in the trace element approximation. 
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1. Introduction 

Since photometric variations were detected in the white 
dwarf HL Tau 76 (Landolt 1968), astronomers have been 
observing multimode pulsations in an increasing num- 
ber of these objects. Of particular interest are the vari- 
able white dwarfs characterized by hydrogen-rich atmo- 
spheres. These variable stars, known in the literature as 
ZZ Ceti or DAV stars, constitute the most numerous 
group amongst degenerate pulsators. Other class of pul- 
sating white dwarfs are the DBV, with helium-rich at- 
mospheres, and the pre-white dwarfs DOVs and PNNVs, 
which show spectroscopically pronounced carbon, oxygen 
and helium features (for reviews of the topic, see Winget 
1988 and Kepler & Bradley 1995). In particular, ZZ Ceti 
stars are found in a narrow interval of effective temper- 
ature (T off ) ranging from 12500 K £ T cfr £ 10700 K. 
Their brightness variations, which reach up to 0.30 mag- 
nitudes, are interpreted as being caused by spheroidal, 
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non-radial g(gravity)-modes of low degree (£ < 2) and 
low and intermediate overtones k (the number of zeros 
in the radial eigenfunction), with periods (Pk) between 2 
and 20 minutes. Radial modes, although found oversta- 
tes in a number of theoretical studies of pulsating DA 
white dwarfs (see, e.g. Saio et al. 1983), have been dis- 
carded as the cause of variability in such stars. This is so 
because the periods involved are shorter than 10 seconds. 
Observationally these high-frequency signatures have not 
been detected thus far. With regard to the mechanism 
that drives pulsations, the k — 7 mechanism is the tradi- 
tionally accepted one (Dolez & Vauclair 1981 and Winget 
et al. 1982). Nonetheless, Brickhill (1991) proposed the 
convective driving mechanism as being responsible for the 
overstability of g-modes in DAVs (see also Goldreich & 
Wu 1999). Although both mechanisms predict roughly 
the observed blue edge of the instability strip, none of 
them arc capable to yield the red edge, where pulsations 
of DA white dwarfs seemingly cease in a very abrupt way 
(Kanaan 1996). 

A longstanding problem in the study of pulsating DA 
white dwarfs is to find the reason of why only a very re- 
duced number of modes are observed, as compared with 
the richness of modes predicted by theoretical studies. 
Indeed, it has been long suspected that some filtering 
mechanism must be acting quite efficiently. The expla- 
nation commonly proposed is that of "mode trapping" 



2 



Corsico et al.: The adiabatic pulsation properties of full DA white dwarf evolutionary models 



phenomenon [j] (Winget et al. 1981; Brassard et al. 1992a; 
Bradley 1996). According to this mechanism, those modes 
(the trapped ones) having a local radial wavelength com- 
parable with the thickness of the hydrogen envelope, re- 
quire low kinetic energy to reach observable amplitudes. 
Then, most of the observed periods should correspond to 
trapped modes. Nevertheless, in a recent asteroseismolog- 
ical study of the ZZ Ceti star G117-B15A (Bradley 1998), 
the observed period of 215 s, which has the larger ampli- 
tude in the power spectrum, does not correspond to the 
trapped mode predicted by the best fitting model. 

The exploration of these very important aspects re- 
quires the construction of detailed DA white dwarf evo- 
lutionary models, particularly regarding the treatment of 
the chemical abundance distribution. Work in this direc- 
tion has recently begun to be undertaken. In fact, Althaus 
et al. (2002) have carried out full evolutionary calculations 
which take into account time dependent element diffusion, 
nuclear burning and the history of the white dwarf pro- 
genitor in a self-consistent way. Specifically, these authors 
have followed the evolution of an initially 3 M© stellar 
model all the way from the stages of hydrogen and helium 
burning in the core through the thermally pulsing and 
mass loss phases till the white dwarf state. Althaus et al. 
(2002) find that the shape of the Ledoux term (an impor- 
tant ingredient in the computation of the Brunt- Vaisala 
frequency; see Brassard et al. 1991, 1992ab) is markedly 
different from that found in previous detailed studies of 
white dwarf pulsations. This is due partly to the effect of 
smoothness in the chemical abundance distribution caused 
by element diffusion, which gives rise to less pronounced 
peaks in the Ledoux term. This, in turn, leads to a sub- 
stantially weaker mode trapping effect, as it has recently 
been found by Corsico et al. (2001). These authors have 
presented the first results regarding the trapping proper- 
ties of the Althaus et al. (2002) evolutionary models. 

The present work is designed to explore at some length 
the pulsation properties of the Althaus et al. (2002) evolu- 
tionary models and to compare our predictions with those 
of others investigators. In addition, the work is intended 
to bring some more insight to the phenomenon of mode 
trapping in the frame of these new evolutionary mod- 
els. In particular, we shall restrict ourselves to analyse 
the same stellar model as that studied in Corsico et al. 
(2001). Specifically, the model, which belongs to the ZZ 
Ceti instability strip, is analyzed in the frame of linear, 
non-radial stellar pulsations in the adiabatic approxima- 
tion. Emphasis will be placed on assessing the role played 
by the internal chemical stratification in the behaviour of 
eigenmodes, and the expectations for the full spectrum 
of periods. Specifically, we shall explore the effects of the 
chemical interfaces on the kinetic energy distribution of 
the modes and their ability to modify the properties of the 
g-mode propagation throughout the star interior. We want 

1 Other physically plausible mechanism of mode filtering is 
related to the interaction of convection with pulsation (see 
Gautschy et al. 1996). 



to mention that, because of the high computational de- 
mands involved in the evolutionary calculations in which 
white dwarf cooling is assessed in a self-consistent way 
with element diffusion and the history of the pre-white 
dwarf, we are forced to restrict our attention to only one 
value for the stellar mass. 

The article is organized as follows. In Sect. ^, we briefly 
describe our evolutionary and pulsational codes. In this 
section we also discuss some aspects concerning the evo- 
lutionary properties of our models. In Sect. we present 
in detail the pulsation results. Finally, Sect. ^| is devoted 
to summarizing our findings. 

2. Details of computations 

Our selected ZZ Ceti model on which the pulsational re- 
sults are based has been calculated by means of an evo- 
lutionary code developed by us at La Plata Observatory. 
The code, which is based on a detailed and up-to-date 
physical description, has enabled us to compute the white 
dwarf evolution in a self-consistent way with the predic- 
tions of time dependent clement diffusion, nuclear burning 
and the history of the white dwarf progenitor. The con- 
stitutive physics includes: new OPAL radiative opacities 
for different metallicities, conductive opacities, neutrino 
emission rates and a detailed equation of state. In addi- 
tion, a network of 30 thermonuclear reaction rates for hy- 
drogen burning (proton-proton chain and CNO bi-cycle) 
and helium burning has been considered. Nuclear reac- 
tion rates are taken from Caughlan & Fowler (1988) and 
Angulo et al. (1999) for the 12 C(a,7) 16 reaction rate. 
This rate is about twice as large as that of Caughlan & 
Fowler (1988). Abundance changes resulting from nuclear 
burning are computed by means of a standard implicit 
method of integration. In particular, we follow the evolu- 
tion of the chemical species 1 H, 3 He, 4 He, 7 Li, 7 Be, 12 C, 
13 C, 14 N, 15 N, 16 0, 17 0, 18 and 19 F. Convection has 
been treated following the standard mixing length theory 
(Bohm-Vitense 1958) with the mixing-length to pressure 
scale height parameter of a =1.5. The Schwarzschild cri- 
tcrium was used to determine the boundaries of convective 
regions. Overshooting and semi-convection were not con- 
sidered. Finally, the various processes relevant for element 
diffusion have also been taken into account. Specifically, 
we considered the gravitational settling, and the chemical 
and thermal diffusion of nuclear species 1 H, 3 He, 4 He, 12 C, 
14 N and 16 O. Element diffusion is based on the treatment 
for multicomponent gases developed by Burgers (1969). 
It is important to note that by using this treatment of 
diffusion we are avoiding the widely used trace element 
approximation (see Tassoul et al. 1990). After computing 
the change of abundances by effect of diffusion, they are 
evolved according to the requirements of nuclear reactions 
and convective mixing. Radiative opacities are calculated 
for metallicities consistent with the diffusion predictions. 
This is done during the white dwarf regime in which grav- 
itational settling leads to metal-depleted outer layers. In 
particular, the metallicity is taken as two times the abun- 
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dance of CNO elements. For more details about this and 
other computational details we refer the reader to Althaus 
et al. (2002) and Althaus et al. (2001). 

We started the evolutionary calculations from a 3 
M Q stellar model at the zero-age main sequence. The 
adopted initial mctallicity is Z= 0.02 and the initial abun- 
dance by mass of hydrogen and helium are, respectively, 
Xn= 0.705 and Xu e = 0.275. Evolution has been com- 
puted at constant stellar mass all the way from the stages 
of hydrogen and helium burning in the core up to the tip of 
the asymptotic giant branch where helium thermal pulses 
occur. After experiencing 11 thermal pulses, the model is 
forced to evolve towards the white dwarf state by invoking 
strong mass loss episodes. The adopted mass loss rate was 
« 10~ 4 M Q yr _1 and it was applied to each stellar model 
as evolution proceeded. After the convergence of each new 
stellar model, the total stellar mass is reduced according 
to the time step used and the mesh points are appropri- 
ately adjusted. As a result of mass loss episodes, a white 
dwarf remnant of 0.563 M Q is obtained. The evolution 
of this remnant is pursued through the stage of planetary 
nebulae nucleus till the domain of the ZZ Ceti stars on 
the white dwarf cooling branch. 

As well known, the shape of the composition transition 
zones plays an important role in the pulsational properties 
of DAV white dwarfs. In this sense, an important aspect 
of these calculations concerns the evolution of the chem- 
ical abundance during the white dwarf regime. In partic- 
ular, element diffusion makes near discontinuities in the 
chemical profile at the start of the cooling branch be con- 
siderably smoothed out by the time the ZZ Ceti domain 
is reached (see Althaus et al. 2002). The chemical pro- 
file throughout the interior of our selected white dwarf 
model is depicted in the upper panel of Fig. 1. Only the 
most abundant isotopes are shown. In particular, the inner 
carbon-oxygen core emerges from the convective helium 
core burning and from the subsequent stages in which the 
helium-burning shell propagates outwards. Note also the 
flat profile of the carbon and oxygen distribution towards 
the centre. This is a result of the chemical rehomogeniza- 
tion of the innermost zone of the star due to Raylcigh- 
Taylor instability (see Althaus et al. 2002). Above the 
carbon-oxygen interior there is a shell rich in both car- 
bon (« 35%) and helium (« 60%), and an overlying layer 
consisting of nearly pure helium of mass 0.003 M©. The 
presence of carbon in the helium-rich region below the 
pure helium layer is a result of the short-lived convective 
mixing which has driven the carbon-rich zone upwards 
during the peak of the last helium pulse on the asymp- 
totic giant branch. We want to mention that the total 
helium content within the star once helium shell burning 
is eventually extinguished amounts to 0.014 M Q and that 
the mass of hydrogen that is left at the start of the cool- 
ing branch is about 1.5 x 10~ 4 M Q , which is reduced to 
7 x 10~ 5 M Q due to the interplay of residual nuclear burn- 
ing and element diffusion by the time the ZZ Ceti domain 
is reached. Finally, we note that the inner carbon-oxygen 
profile of our models is somewhat different from that of 



Salaris et al. (1997). In particular, we find that the size 
of the carbon-oxygen core is smaller than that found by 
Salaris et al. (1997) with the consequent result that the 
drop in the oxygen abundance above the core is not so 
pronounced as in the case found by these authors. We 
suspect that this different behaviour could be a result of a 
different treatment of the convective boundary during the 
core helium burning. 

For the pulsation analysis we have employed the code 
described in Corsico & Benvenuto (2002). We refer the 
reader to that paper for details. Here we shall describe 
briefly our strategy of calculation and mention the pulsa- 
tion quantities computed that are relevant in this study. 
The pulsational code is based on the general Newton- 
Raphson technique (like the Henyey method employed in 
stellar evolution studies). The code solves the differential 
equations governing the linear, non-radial stellar pulsa- 
tions in the adiabatic approximation (see Unno et al. 1989 
for details of their derivation). The boundary conditions at 
the stellar centre and surface are those given by Osaki & 
Hansen (1973) (see Unno et al. 1989 for details). Following 
previous studies of white dwarf pulsations, the normaliza- 
tion condition adopted is Sr/r = 1 at the stellar surface. 
After selecting a starting stellar model we choose a con- 
venient period window and the interval of interest in T e g. 
The evolutionary code computes the white dwarf cooling 
until the hot edge of the T e g-interval is reached. Then, the 
program calls the set of pulsation routines to begin the 
scan for modes. In order to obtain the first approxima- 
tion to the eigenfunctions and the eigenvalue of a mode 
we have applied the method of the discriminant (Unno 
et al. 1989). Specifically, we adopt the potential bound- 
ary condition (at the surface) as the discriminant function 
(see Corsico & Benvenuto 2002). When a mode is found, 
the code generates an approximate solution which is it- 
eratively improved to convergence (of the eigenvalue and 
the eigenfunctions simultaneously) and then stored. This 
procedure is repeated until the period interval is covered. 
Then, the evolutionary code generates the next stellar 
model and calls pulsation routines again. Now, the pre- 
viously stored modes are taken as initial approximation 
to the modes of the present stellar model and iterated to 
convergence. 

For each computed mode we obtain the cigenperiod 
Pk (Pk — 2^/(7^, being <7k the eigenfrequency) and the di- 
mcnsionless eigenfunctions j/i, • ■ ■ , J/4 (see Unno et al. 1989 
for their definition). With these eigenfunctions and the di- 
mensionless eigenvalue lo\ — <r^(GM*/ R^) -1 we compute 
for each mode considered the oscillation kinetic energy, 
(-Ekin)k, given by: 



(£kin)fc - ^{GM*R\ 



in f 

Jo 



x 2 p 



2 2 , 2 1(1+1) 2 
X wl + X In . .2^,2^2 



dx,(l) 



and the first order rotation splitting coefficient, C^k, 
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where M* and P* are the stellar mass and the stellar 
radius respectively, G is the gravitation constant, C\ = 
(r I -R*) 3 (M» / M r ) and x = r/P*. In addition, we compute 
the weight functions, WF, and the variational period, P£ , 
as given by Kawaler et al. (1985). Finally, for each model 
computed we derive the asymptotic spacing of periods, 
AP A , given by (Tassoul 1980; Tassoul et al. 1990): 



AP A = 



and Pq is denned as 



Pn = 2tt z 



X2 N 

— ax 

x 



(3) 



(4) 



where X2 correspond to the location of the base of the 
outer convection zone. The Brunt- Vaisala frequency (TV), 
a fundamental quantity of white dwarf pulsations, is com- 
puted employing the "modified Ledoux" treatment. This 
treatment explicitly accounts for the contribution to N 2 
from any change in composition in the interior of model 
(the zones of chemical transition) by means of the Ledoux 
term B (see Brassard et al. 1991). We want to mention 
that we have also employed a numerical differentiation 
scheme for computing N 2 directly from its definition. We 
found that this scheme yields the same results as those 
derived from the modified Ledoux treatment. 

As mentioned, for the pulsation analysis in this study 
we have selected a white dwarf model representative of the 
ZZ Ceti instability band. Specifically, we have picked out 
a 0.563 M model at T cff w 12000 K. In Table 1 we show 
the main characteristics of our template model. In the in- 
terests of comparison, the same quantities corresponding 
to a similar model of P. Bradley (2002) (private communi- 
cation) are shown. The Brunt- Vaisala frequency and the 
Ledoux term (P) corresponding to our template model are 
shown in the middle and bottom panels of Fig. 1 in terms 
of the outer mass fraction. Note the particular shape of 
P, which is a 



Table 1. 




Our template 


Bradley's model 




model 




M*/M 


0.563 


0.560 


Tcff [K] 


11996 


12050 


log(L*/L ) 


-2.458 


-2.462 


log(P»/R ) 


-1.864 


-1.866 


log(M H /M») 


-3.905 


-3.824 


log(M He /M*) 


-1.604 


-1.824 


k)g(p c ) [g cm" 3 ] 


6.469 


6.466 


log(P c ) [10 6 K] 


7.086 


7.087 



direct consequence of the chemical profile. In turn, the 
features of B are reflected in the Brunt- Vaisala frequency. 
As a result of element diffusion, the chemical profiles of 
our evolutionary models are very smooth in the inter- 
faces (see upper panel of Fig. 1), and this explains the 
presence of extended tails in the shape of B. Also, note 



that our model is characterized by a chemical interface in 
which three ionic species in appreciable abundances co- 
exists: helium, carbon and oxygen (see early in this sec- 
tion). This transition region gives two contributions to 
B, one of them of relatively great magnitude, placed at 
logg ~ —1.4, and the other, more external and of very 
low height at \ogq ~ —2.2. As a last remark, the contri- 
bution of the hydrogen-helium interface to B is less than 
that corresponding to the helium-carbon-oxygen transi- 
tion. Note that the contributions from the Ledoux term 
are translated into smooth bumps on N 2 . The characteris- 
tics of B and iV 2 as predicted by our models are markedly 
different from those found in previous studies in which the 
white dwarf evolution is treated in a simplified way, par- 
ticularly regarding the chemical abundance distribution in 
the outer layers (e.g., Tassoul et al. 1990; Brassard et al. 
1991, 1992ab; Bradley 1996). For more details see Althaus 
et al. (2002). 



3. Results 

For our template model we have computed g-modes with 
£=1,2 and 3 (we do this because geometric cancellation 
effects grow progressively for larger I in non-radial oscil- 
lations; see Dziembowski 1977), with periods in the range 
of 50 s <; Pfc <; 1300 s. Let us quote that for mode cal- 
culations we have employed up to 5000 mesh points. For 
all of our pulsation calculations, the relative difference be- 
tween Pfc and P)! remains lower than 10~ 3 . This gives an 
indication of the accuracy of our calculations. 

We begin by examining Figs. 2 to 4, the upper panels 
of which show the logarithm of the oscillation kinetic en- 
ergy of modes with, respectively, £ = 1,2 and 3 in terms 
of computed periods. Middle panels depict the values for 
the forward period spacing AP& (= Pfc+i — P&) together 
with the asymptotic value APa as given by dotted lines 
0. Finally, in the bottom panel of these figures we de- 
pict the Cijt values as well as the asymptotic values (dot- 
ted lines) that these coefficients adopt for high overtones, 
that is Cg )k w !/£(£ + 1) (Brickhill 1975). An inspection 
of plots reveals some interesting characteristics. To begin 
with, the quantities plotted exhibit two clearly different 
trends. Indeed, for Pf. > 500 — 600 s and irrespective 
of the value of £, the distribution of oscillation kinetic 
energy is quite smooth. Note that the log(Pkin)fc values 
of adjacent modes are quite similar, which is in contrast 
with the situation found for lower periods. On the other 
hand, the period spacing diagrams show appreciable de- 
partures of APfc from the asymptotic prediction (Eq. 3) 
for Pfc < 500 — 600 s. As well known, this is due mostly 
to the presence of chemical abundance transitions in DA 
white dwarfs. In contrast, for higher periods the APfc of 



2 Strictly speaking, APa as given by Eq. Q corresponds to 
the asymptotic (high k) separation of periods for chemically 
homogeneous and radiative stars. However, it is very close to 
the mean period spacing even for chemically stratified white 
dwarfs. See Tassoul et al. (1990) for details. 
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the modes tend to APa- Also, note that the Cg^ values 
tend to the asymptotic value for Pk ^ 500 — 600 s. 

An important aspect of the present study is related to 
the mode trapping and confinement properties of our mod- 
els. For the present analysis we shall employ the weight 
functions, WF. We elect WF because this function gives 
the relative contribution of the different regions in the star 
to the period formation (Kawaler et al. 1985; Brassard et 
al. 1992ab). We want to mention that we have also care- 
fully examined the density of kinetic energy (the integrand 
of Eq. 1) for each computed mode. For our purposes here, 
this quantity gives us basically the same information that 
provided by WF. We show in Fig. 5 to 7 the WF for all 
of the computed modes corresponding to I = 1. In addi- 
tion, we include in each plot of these figures the Ledoux 
term in arbitrary units (dotted lines) in order to make 
easier the location of the chemical transition regions of 
the model. In the interests of a proper interpretation of 
these figures, we suggest the reader to see also Fig. 2. For 
low periods, a variety of behaviour is encountered. For 
instance, the gi mode is characterized by a WF corre- 
sponding to the well known mode trapping phenomenon, 
that is, gi is formed in the very outer layers irrespective 
of the details of the deeper chemical profile, as previously 
reported by previous studies (see Brassard et al. 1992ab). 
In contrast, it is the helium-carbon-oxygen transition that 
mostly contributes to the formation of the g2 mode, whilst 
the hydrogen-helium transition plays a minor role. This 
mode would be representative of the "confined modes" ac- 
cording to Brassard et al. (1992ab). WF for modes g3 and 
g4 is qualitatively similar to that of gi, except that they 
are not exclusively formed in the hydrogen-rich envelope, 
but also in the helium-carbon-oxygen interface. On the 
other hand, the high-density zone underlying the helium- 
carbon-oxygen transition plays a major role in the forma- 
tion of mode gs. From Eq. (1) is clear that the (-Ekin)fc V& 1~ 
ues are proportional to the integral of the squared eigen- 
functions, weighted by p. As a result, the gs mode is char- 
acterized by a high oscillation kinetic energy value (see 
Fig. 2). Note that the helium-carbon-oxygcn transition 
region also contributes to the formation of modes g6 and 
gio- The gio mode is particularly interesting, because it is 
formed over a wide range of the stellar interior, thus being 
also a high kinetic energy mode. The WFs correspond- 
ing to remaining modes do not differ appreciably amongst 
them. They exhibit contributions mainly from the outer 
layers of the model, though they also show small ampli- 
tudes in deeper regions. Note that for all of the modes 
shown in Figs. 5 to 7 there is a strong contribution to 
WFs from the hydrogen-helium transition region. This 
indicate that, as found in previous studies, this chemical 
interface plays a fundamental role in the period formation 
of modes. We want to mention that we have elected for 
this analysis the dipolar modes (£ = 1) for brevity; the 
results for I = 2,3 arc qualitatively similar to those of 
1=1. 

From the analysis performed above based on the 
weight functions, we can clearly appreciate that for Pk <; 



500 — 600 s the outer layer of the model appreciably con- 
tribute to the WFs. This is expected, because, as well 
known, (/-modes in white dwarfs are envelope modes. As 
mentioned, the WFs of high order modes are very simi- 
lar, indicating that these modes have essentially the same 
characteristics. At this point, we could, in principle, clas- 
sify these modes either like trapped or partially trapped in 
the outer envelope or like "normal" modes (in the termi- 
nology of Brassard et al. 1992), that is, without enhanced 
or diminished oscillation kinetic energies as in the case 
of eigenmodes corresponding to chemically homogeneous 
stellar models. In fact, the curves log^kin)*: — Pk depicted 
in Figs. 2 to 4, in particular for periods exceeding 500 — 600 
s, strongly resemble the kinetic energy distribution cor- 
responding to a model in which there no exist chemical 
interfaces. With the aim of solving such an ambiguity, we 
have performed pulsation calculations arbitrarily setting 
B = in the computation of the Brunt- Vaisala frequency. 
As mentioned, the modified Ledoux treatment employed 
in the computation of the Brunt- Vaisala frequency bears 
explicitly the effect from changes in chemical composition 
by means of the Ledoux term B. So, by forcing B = 
the effects of the chemical transitions are strongly min- 
imized (but not completely eliminated; see inset of Fig. 
lc) on the whole pulsational pattern. In this way we ob- 
tain an approximate chemically homogeneous white dwarf 
model (see Brassard et al. 1992b for a similar numerical 
experiment). The oscillation kinetic energy values result- 
ing for this simulated "homogeneous" model are shown 
in Fig. 8 with dotted lines. In the interests of a com- 
parison, we show the results corresponding to our (full) 
template model with solid lines. It is clear from the fig- 
ure that the distribution of (-Ekin)fc values for both sets 
of computations (and for each value of t) is very simi- 
lar in the region of long periods. However, note that the 
curves corresponding to the modified model are shifted to 
higher energies (by « 0.2 dcx) as compared with the sit- 
uation of the full model. We have carefully compared the 
WF for each mode of the full model with the correspond- 
ing mode of the "homogeneous" model (i.e. modes which 
have closest period values although generally for different 
radial order k). We found that, for modes with periods 
exceeding « 600 s, the WFs are almost identical in both 
cases at the regions above the hydrogen-helium transition. 
However, below this interface the WFs corresponding to 
the "homogeneous" model show larger amplitudes as com- 
pared with the case of the full model. Thus, we can con- 
clude that for the full model, all the modes correspond- 
ing to the long period region of the pulsational spectrum 
must be considered as partially trapped in the hydrogen- 
rich envelope. In others words, the chemical distribution 
at the hydrogen-helium transition has noticeable effect on 
each mode, but this effect is the same for all modes. This 
conclusion is reinforced by the fact that the first order ro- 
tation splitting coefficients (Ci.k) for the full model adopt 
higher values as compared with those corresponding to the 
"homogeneous" model (not shown here for brevity), thus 
lying nearest to the asymptotic prediction. As found by 
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Brassard et al. (1992ab), it is an additional characteris- 
tic feature of trapped modes in the hydrogen-rich outer 
region of white dwarfs. 

An important finding of this work is the effect of chemi- 
cal abundance distribution resulting from time dependent 
diffusion on the mode trapping properties in DA white 
dwarfs. In fact, as shown in Fig. 2 to 4, for periods exceed- 
ing « 500 — 600 s, the distribution of (.Ekin)fc is smooth, 
and APfc values tend to the asymptotic value. This is quite 
different from that found in previous studies. Our calcu- 
lations reveal that the capability of mode filtering due 
to mode trapping effects virtually vanish for high peri- 
ods when account is made of white dwarf models with 
diffusively evolving chemical stratifications (see Corsico 
et al. 2001). In order to make a detailed comparison of 
the predictions of our models with those found in previ- 
ous studies we have carried out additional pulsational cal- 
culations by assuming diffusive equilibrium in the trace 
element approximation at the hydrogen-helium interface 
(see Tassoul et al. 1990). This treatment has been com- 
monly invoked in most of the pulsation studies to model 
the composition transition regions. The resulting hydro- 
gen chemical profile and the corresponding Ledoux term 
and Brunt- Vaisala frequency TV are shown in Fig. 9, to- 
gether with the predictions of time dependent element dif- 
fusion. The trace element assumption leads to an abrupt 
change in the slope of the chemical profile which is re- 
sponsible for the pronounced peak in the Brunt- Vaisala 
frequency at log(l — M r /M») w —4. As can be clearly 
seen in Fig. 10 for £ = 1 to 3, the diffusive equilibrium in 
the trace element approximation gives rise to an oscillation 
kinetic energy spectrum and period spacing distribution 
that are substantially different from those given by the 
full treatment of diffusion (see Figs. 2 to 4), particularly 
for high periods. The most outstanding feature depicted 
by Fig. 10 is the trapping signatures exhibited by certain 
modes both in the log(Pkin)fc and AP& values. This is in 
agreement with other previous results (see Brassard et al. 
1992b, particularly their figures 20a and 21a for the case 
of Mh = lO^Af*)^. As well known, trapped modes corre- 
spond to those modes which are characterized by minima 
in their oscillation kinetic energy values and local minima 
in the period spacing having the same fc-value or differing 
by 1. For the purpose of illustration, we compare in Fig. 11 
and 12 the predictions of equilibrium diffusion in the trace 
element approximation and time-dependent element diffu- 
sion, respectively, for WF corresponding to the modes g38 
and g39 with £ = 2. Clearly, in the case of diffusive equi- 
librium in the trace element approximation, mode g3g cor- 
responds to a trapped mode characterized by small values 



3 Note that, as found by Brassard et al. (1992a), the contrast 
between the kinetic energies values of the trapped modes and 
those of the non-trapped ones is not very large for massive hy- 
drogen envelopes (like that presented in this work), due to the 
fact that the hydrogen- helium interface is located in a deep, 
highly degenerate region, where the amplitudes of eigenfunc- 
tions are very small. 



of the weight function below the hydrogen-helium transi- 
tion, as compared with the adjacent, non-trapped mode 
g 38 . By contrast, such modes show very similar amplitudes 
of their WF when account is taken of a full diffusion treat- 
ment to model the composition transition regions (see Fig. 
12). We would also like to comment on the fact that the 
diffusive equilibrium condition is far from being reached 
at the bottom of the hydrogen envelope of our model. In 
Althaus et al. (2002) we argued that the situation of diffu- 
sive equilibrium in the deep layers of a DAV white dwarf 
is not an appropriate one for describing the shape of the 
chemical composition at the hydrogen-helium transition 
zone. In fact, during the ZZ Ceti stage time-dependent 
diffusion modifies the spatial distribution of the elements, 
particularly at the chemical interfaces (see also Iben & 
MacDonald 1985). In addition, for the case of thick hy- 
drogen envelopes, we have recently found that under the 
assumption of diffusive equilibrium, a white dwarf does 
not evolve along the cooling branch, but rather it expe- 
riences a hydrogen thermonuclear shell flash (see Corsico 
et al. 2002). This is so because if diffusion had plenty of 
time to evolve to an equilibrium situation then the tail of 
the hydrogen distribution would have been able to reach 
hot enough layers to be ignited in a flash fashion. 

To place some of the results of the foregoing paragraph 
in a more quantitative basis, we list in Tables 2 and 3 the 
values for P&, APfc and log(Pkin)fc for modes correspond- 
ing to £ = 1 and £ = 2, in the case of equilibrium dif- 
fusion and time dependent element diffusion. In Table 2, 
the "m" corresponds to minima, and "M" stands for max- 
ima. We have labeled the minima of AP& and the minima 
and maxima of log(Pkin)fc, in correspondence with Fig. 10. 
Note that for the case with £ = 1 there is a direct correla- 
tion (indicated by arrows) between minima in log(Pkin)fc 
and APfc_i for most of high order modes, whereas for the 
case with £ = 2 this correspondence is between minima in 
log(-Ekin)fe and Apj. The modes with minima in kinetic 
energy are classified as trapped (T) ones. In contrast to 
the case of equilibrium diffusion, the results correspond- 
ing to the time dependent element diffusion treatment do 
not show clear minima or maxima in kinetic energy, as 
can be appreciated in Figs. 2 to 4 and Table 3. We have 
compared the periods of our model with those kindly pro- 
vided by Bradley corresponding to his 0.560 M & white 
dwarf model, and we find that our periods are typically 6 
% shorter. In part, this difference is due to the somewhat 
smaller mass of the Bradley's model and the different in- 
put physics characterizing both stellar models. 



On the basis of our results, we claim that the treatment 
of the chemical profile at the chemical transitions is a key 
ingredient in the computation of the g-spectrum. This is 
particularly true regarding the mode trapping properties, 
which are considerably altered when a physically sound 
treatment of the chemical evolution is incorporated in such 
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Table 2. Pulsation properties for £ = 1, 2 modes 
corresponding to the case of equilibrium diffusion. 



1 = 1 


£ = 2 






AP A = 46.26 


s 








AP A = 26.71 


s 




k 


Pk 




AP k 




log ( Skin )fe 


k 


Pk 




AP k 




log(_Ekin)fe 




[s] 




L J 




[erg] 




[s] 




L J 




[erg] 


1 


126.99 


T 


76.38 




45.84 m 


1 


73.39 


T 


48.02 




45.84 m 


2 


203.37 




73.91 




46.97 M 


2 


121.40 




39.77 




46.85 M 


3 


277.28 


T 


26.02 m 




44.60 m 


3 


161.17 


T 


18.11 m 




44.52 m 


4 


303.30 




40.79 




44.61 m 


4 


179.28 




39.51 




44.61 M 


5 


344.09 




40.65 




44.93 M 


5 


218.80 


T 


14.17 m 




43.77 m 


6 


384.75 




44.00 




43.62 m 


6 


232.97 




17.06 




43.89 M 


7 


428.75 




61.25 




43.43 


7 


250.03 




33.50 




43.49 


8 


490.00 


T 


59.93 




42.90 m 


8 


283.53 




40.06 




42.89 


9 


549.93 




21.12 m 


/ 


42.99 M 


9 


323.59 


T 


24.67 




42.53 m 


10 


571.05 




47.22 




42.59 


10 


348.26 




18.24 m 


/ 


42.48 M 


11 


618.27 




44.54 m 




42.23 


11 


366.50 




19.88 




42.27 


12 


662.81 


T 


48.65 


\ 


41.91 m 


12 


386.39 




35.55 




41.95 


13 


711.46 




30.69 m 




42.24 M 


13 


421.94 


T 


21.12 m 


( 


41.76 m 


14 


742.15 


T 


46.66 


\ 


41.75 m 


14 


443.06 




30.37 




41.96 M 


15 


788.81 




42.28 m 




41.98 M 


15 


473.43 


T 


19.88 m 




41.58 m 


16 


831.09 


T 


48.51 


\ 


41.48 m 


16 


493.32 




25.75 




41.66 M 


17 


879.60 




43.99 m 




41.55 M 


17 


519.06 




19.59 m 




41.45 


18 


923.59 


T 


48.15 


\ 


41.16 m 


18 


538.65 




30.84 




41.21 


19 


971.75 




33.31 m 




41.31 M 


19 


569.49 


T 


22.68 m 




41.10 m 


20 


1005.05 


T 


47.02 


\ 


41.16 m 


20 


592.18 




27.64 




41.29 M 


21 


1052.07 




42.38 m 




41.32 M 


21 


619.82 




19.68 m 




41.22 


22 


1094.45 


T 


55.95 


\ 


41.10 m 


22 


639.51 




32.53 




41.26 M 


23 


1150.40 




32.44 m 




41.22 M 


23 


672.04 


T 


24.88 m 




41.00 m 


24 


1182.83 


T 


43.15 


\ 


41.09 m 


24 


696.92 




27.97 




41.12 M 


25 


1225.98 




40.36 m 




41.16 M 


25 


724.89 


T 


21.88 m 




40.93 m 


26 


1266.34 


T 


53.93 


\ 


40.93 m 


26 


746.77 




25.83 




41.25 M 


27 


1320.27 




38.87 m 




41.11 M 


27 


772.60 


T 


23.68 m 




41.01 m 














28 


796.28 




28.33 




41.18 M 














29 


824.61 


T 


22.94 m 




41.03 m 














30 


847.55 




27.53 




41.17 M 














31 


875.08 


T 


21.89 m 




41.04 m 














32 


896.97 




30.30 




41.14 M 














33 


927.27 


T 


25.47 m 


^_ 


40.98 m 














34 


952.74 




27.78 




41.22 M 














35 


980.52 


T 


23.58 m 


( 


41.06 m 














36 


1004.10 




26.44 




41.32 M 














37 


1030.54 


T 


23.93 m 




41.12 m 














38 


1054.46 




28.40 




41.35 M 














39 


1082.87 


T 


23.20 m 




41.18 m 














40 


1106.07 




26.52 




41.48 M 














41 


1132.58 


T 


23.99 m 


<_ 


41.27 m 














42 


1156.57 




29.72 




41.46 M 














43 


1186.29 


T 


25.32 m 




41.32 m 














44 


1211.61 




27.85 




41.57 M 














45 


1239.47 


T 


24.74 m 




41.41 m 














46 


1264.20 




26.62 




41.68 M 














47 


1290.82 


T 


25.45 m 




41.48 m 














48 


1316.27 




26.83 




41.70 M 



calculations. This conclusion is valid at least for massive 
hydrogen envelopes as predicted by our full evolutionary 
calculations. To get a deeper insight into these aspects, we 
examine the node distribution of the eigenfunctions (Figs. 
13 and 14). According to Brassard et al. (1992a), this 
is an useful diagnostic for mode trapping. Here, a mode 



is trapped above the hydrogen-helium interface when its 
eigenfunction y\ has a node just above of such an inter- 
face, and the corresponding node in yi lies just below that 
interface. Note that this statement is clearly satisfied by 
our model with diffusive equilibrium in the trace element 
approximation, as shown in Fig. 13. In this figure, the ver- 
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Table 3. Same as Table 2, but for the case of time-dependent 
element diffusion. 



1 = 1 


1 = 2 




AP A 


= 45.83 


s 




APa 


= 26.30 s 


k 


Pk 


AP fc 


log(P kin ) fc 


k 


Pk 


AP k 


log(-Ekin)fc 




[s] 


[s] 


[erg] 




[s] 


[s] 


[erg] 


1 


126.98 


76.87 


45.84 


1 


72.14 


49.34 


45.84 


2 


203.85 


66.57 


46.83 


2 


121.48 


34.89 


46.69 


3 


270.43 


33.77 


44.77 


3 


156.38 


21.63 


44.72 


4 


304.19 


35.98 


44.48 


4 


178.00 


32.08 


44.41 


5 


340.17 


36.30 


44.69 


5 


210.09 


23.30 


44.10 


6 


376.47 


44.67 


43.97 


6 


233.39 


10.74 


43.84 


7 


421.13 


67.75 


43.30 


7 


244.12 


36.72 


43.40 


8 


488.88 


51.08 


43.04 


8 


280.85 


31.47 


43.05 


9 


539.95 


26.84 


42.78 


9 


312.31 


33.28 


42.67 


10 


566.80 


46.27 


43.04 


10 


345.59 


16.62 


42.44 


11 


613.06 


44.76 


42.18 


11 


362.21 


19.97 


42.30 


12 


657.83 


46.99 


42.03 


12 


382.18 


30.40 


42.05 


13 


704.82 


32.45 


42.08 


13 


412.58 


25.26 


41.87 


14 


737.27 


42.44 


41.93 


14 


437.84 


26.30 


41.86 


15 


779.71 


44.78 


41.88 


15 


464.14 


24.63 


41.78 


16 


824.49 


43.42 


41.68 


16 


488.77 


21.89 


41.58 


17 


867.91 


50.43 


41.50 


17 


510.66 


22.90 


41.59 


18 


918.34 


43.07 


41.23 


18 


533.56 


24.96 


41.30 


19 


961.42 


38.91 


41.20 


19 


558.53 


27.11 


41.14 


20 


1000.33 


40.24 


41.28 


20 


585.64 


24.25 


41.17 


21 


1040.57 


45.71 


41.19 


21 


609.89 


22.92 


41.21 


22 


1086.28 


50.32 


41.19 


22 


632.80 


26.05 


41.27 


23 


1136.60 


40.88 


41.13 


23 


658.86 


28.13 


41.09 


24 


1177.48 


36.91 


41.17 


24 


686.98 


27.00 


41.00 


25 


1214.40 


42.99 


41.06 


25 


713.98 


23.95 


40.97 


26 


1257.38 


46.02 


40.96 


26 


737.94 


22.71 


41.01 


27 


1303.40 


46.39 


40.99 


27 


760.64 


24.69 


41.07 










28 


785.33 


25.47 


41.02 










29 


810.80 


25.67 


41.05 










30 


836.47 


24.62 


41.04 










31 


861.09 


24.42 


41.04 










32 


885.51 


25.09 


41.05 










33 


910.60 


27.09 


41.02 










34 


937.69 


27.41 


41.03 










35 


965.11 


25.02 


41.09 










36 


990.12 


24.63 


41.15 










37 


1014.75 


24.66 


41.18 










38 


1039.41 


25.64 


41.21 










39 


1065.05 


26.07 


41.23 










40 


1091.12 


24.20 


41.29 










41 


1115.32 


24.33 


41.35 










42 


1139.66 


26.34 


41.36 










43 


1166.00 


27.19 


41.37 










44 


1193.19 


26.48 


41.42 










45 


1219.67 


26.20 


41.46 










46 


1245.87 


25.35 


41.52 










47 


1271.22 


25.13 


41.56 










48 


1296.35 


25.98 


41.58 



tical dotted line at r = 0.927-R* indicate the location of the 
pronounced peak in the Brunt- Vaisala frequency (see Fig. 
9). However, the node distribution becomes markedly dif- 
ferent in the sequence with non-equilibrium diffusion and 
does not seem possible to find a well defined interface in 
that case (see Fig. 14). Thus, it does not seem to be clear 



that the above-mentioned trapping rule can be directly 
applied to this case. In fact, because the hydrogen-helium 
interface becomes very smooth in our models, the peak 
in the Brunt- Vaisala frequency is not very pronounced. 
Accordingly, the capability of mechanical resonance of our 
model turns out to be weaker. This causes the node dis- 
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tribution in the cigenfunctions to be quite different from 
that corresponding to the diffusive equilibrium approach. 



4. Conclusions 

In this work we have explored the pulsational proper- 
ties of detailed evolutionary models recently developed by 
Althaus et al. (2002). Attention has been focused on a 
ZZ Ceti model in the frame of linear, non-radial oscilla- 
tions in the adiabatic approximation. White dwarf cool- 
ing has been computed in a self-consistent way with the 
evolution of the chemical abundances resulting from the 
various diffusion processes and nuclear burning. Element 
diffusion is based on a multicomponent gas treatment; so, 
the trace element approximation is avoided in our calcu- 
lations. In addition, the evolutionary stages prior to the 
white dwarf formation have been considered. In partic- 
ular, element diffusion causes near discontinuities in the 
chemical profile at the start of the cooling branch to be 
considerably smoothed out by the time the ZZ Ceti do- 
main is reached. 

An important aspect of this work has been to assess 
the role played by the internal chemical stratification of 
these new models in the behaviour of the eigenmodes, and 
the expectations for the full g-spectrum of periods. We 
have analyzed mainly the mode weight functions, which 
show the regions of the star that mostly contribute to 
the period formation. Our study suggests the existence of 
a much wider diversity of eigenmodes for periods shorter 
than « 500 — 600 s than found in previous works (Brassard 
et al. 1992ab). An important finding of this study is the 
effect of time-dependent element diffusion on the mode 
trapping properties in DA white dwarfs. We conclude that 
for periods longer than « 500 — 600 s all of the modes 
seem to be partially trapped in the hydrogen-rich enve- 
lope of the star. This conclusion, based on the fact that 
the weight functions of these modes show low amplitudes 
(WF <J 0.1) below the hydrogen-helium transition (even 
lower as compared with the case of a simulated chemically 
homogeneous model), implies that the capability of mode 
selection due to mode trapping effects vanishes for high 
periods when account is made of white dwarf models with 
diffusively evolving stratifications. This conclusion is valid 
at least for massive hydrogen envelopes as predicted by our 
full evolutionary calculations. This behaviour is markedly 
different from that found in other studies based on the 
assumption of diffusive equilibrium in the trace element 
approximation. This assumption leads to a pronounced 
peak in the Ledoux term at the hydrogen-helium inter- 
face, which is responsible for the trapping of modes in the 
outer hydrogen-rich layers. We have verified this fact by 
performing additional pulsation calculations on a model in 
which the chemical profile at the hydrogen-helium transi- 
tion is given by equilibrium diffusion in the trace element 
approximation. Finally, the prediction of both diffusion 
treatments for the node distribution of the eigenfunctions 
has been compared. We found that node distribution at 



the hydrogen-helium chemical interface is very sensitive 
to the treatment of the chemical profile at that interface. 

On the basis of these new results, we are forced to con- 
clude that for high periods, trapping mechanism in mas- 
sive hydrogen envelopes of stratified DA white dwarfs is 
not an appropriate one to explain the fact that all the 
modes expected from theoretical models are not observed 
in ZZ Ceti stars. Interestingly, a weaker trapping effect 
on the periodicities in DB white dwarfs has also been re- 
ported (Gautschy & Althaus 2001). We think that the 
results presented in this work deserves further exploration 
from the point of view of a non-adiabatic stability analy- 
sis. Work in this direction is in progress. 
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Fig. 1. Panel a: The chemical abundance distribution of our stellar model for hydrogen (dotted line), helium (dashed 
line), carbon (solid line) and oxygen (dot-dashed line). Panel b: the Ledoux term, B. Panel c: the logarithm of the 
squared of the Brunt- Vaisala, frequency (^V 2 ). In the inset of this panel, the logarithm of the squared of the Brunt- 
Vaisala, frequency computed neglecting the term B is depicted. Note that the imprints of the chemical transition zones 
in the functional form of N 2 are not completely eliminated, in particular at the helium-carbon-oxygen and hydrogen- 
helium interfaces. The stellar mass of the white dwarf model is 0.563 M© and the effective temperature is w 12000K. 
Quantities are shown in terms of the outer mass fraction. 
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Fig. 2. The logarithm of the oscillation kinetic energy, forward period spacing and first order rotation splitting co- 
efficient (upper, middle and lower panels, respectively) for modes with I = 1 as a function of computed periods. 
Dotted lines show the asymptotic behaviour for APk and Cg y k- The values of (-Eidn)fc correspond to the normalization 
j/i = Sr/r = 1 at r = R*. 
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Fig. 3. Same as Fig. 2, but for 1 = 2. 
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Fig. 4. Same as Fig. 2, but for I = 3. 
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Fig. 5. The normalized weight function (solid lines) in terms of the outer mass fraction, for modes gi to gs with I = 1. 
In the interests of comparison, dotted lines depict the run of the Ledoux term (arbitrary units). 
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Fig. 6. Same as Fig. 5, but for modes gg to gi6- 
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Fig. 7. Same as Fig. 5, but for modes gi7 to g24- 
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Fig. 8. The logarithm of the oscillation kinetic energy for modes with I = 1, £ = 2 and I = 3 (upper, middle and lower 
panel, respectively), as a function of computed periods. In the interests of clarity, symbols corresponding to eigenmodes 
have been omitted. Solid lines correspond to our template model, and dotted lines correspond to the "homogeneous" 
model in which 5 = 0. 
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Fig. 9. Upper panel: hydrogen abundance distribution at the hydrogen-helium interface as given by multicomponent, 
non-equilibrium diffusion (solid line) and diffusive equilibrium in the trace element approximation (thin line). Lower 
panel: the logarithm of the squared Brunt- Vaisala frequency for the both treatment of diffusion mentioned above. The 
inset shows the prediction for the Ledoux term. For details, see text. 
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Fig. 10. The logarithm of the oscillation kinetic energy (upper panel) and period spacing (lower panel) for 1=1 
(figure a), £ = 2 (figure b) and £ = 3 (figure c) in terms of the computed periods, for the case of diffusive equilibrium 
in the trace element approach. As found in previous studies, this approximation gives rise to a kinetic energy pattern 
and spacing of consecutive periods in which trapping signatures are clearly noted. This is in contrast to the prediction 
of time dependent element diffusion given in Figs. 2 to 4. 



Corsico et al.: The adiabatic pulsation properties of full DA white dwarf evolutionary models 



21 




Log(l-M r /M,) 



Fig. 11. The normalized weight function for modes g38 and g3g with I = 2 corresponding to the stellar model in which 
the hydrogen-helium chemical transition has been treated in the diffusive equilibrium and trace element approximation. 
Note the lower amplitude of WF below the hydrogen-helium transition for mode with k— 39, which corresponds to a 
trapped one in the hydrogen envelope. 
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Fig. 12. Same of Fig. 11, but for the case in which the hydrogen-helium chemical transition has been computed 
assuming time dependent element diffusion. 
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Fig. 13. The node distribution of cigenfunctions j/i (filled dots) and (empty dots) at the hydrogen-helium transition 
region for modes with t = 2, according to the diffusive equilibrium prediction in the trace element approximation. 
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Fig. 14. Same as Fig. 13, but for the case of chemical profiles resulting from time dependent element diffusion. See 
text for details. 



